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ABSTRACT 

We  have  recently  discovered  a  supermassive  binary  black  hole  system  with  a  projected  separation  between  the 
two  black  holes  of  7.3  pc  in  the  radio  galaxy  0402+379  (Rodriguez  et  al.  2006).  This  is  the  most  compact 
supermassive  binary  black  hole  pair  yet  imaged  by  more  than  two  orders  of  magnitude.  We  present  Global 
VLBI  observations  at  1.3464  GHz  of  this  radio  galaxy,  taken  to  improve  the  quality  of  the  Hi  data.  Two 
absorption  lines  are  found  toward  the  southern  jet  of  the  source,  one  redshifted  by  370  zb  10  km  s-1  and  the 
other  blueshifted  by  700  ±10  km  s-1  with  respect  to  the  systemic  velocity  of  the  source,  which,  along  with 
the  results  obtained  for  the  opacity  distribution  over  the  source,  suggests  the  presence  of  two  mass  clumps 
rotating  around  the  central  region  of  the  source.  We  propose  a  model  consisting  of  a  geometrically  thick  disk, 
of  which  we  only  see  a  couple  of  clumps,  that  reproduces  the  velocities  measured  from  the  Hi  absorption 
profiles.  These  clumps  rotate  in  circular  Keplerian  orbits  around  an  axis  that  crosses  one  of  the  supermassive 
black  holes  of  the  binary  system  in  0402+379.  We  find  an  upper  limit  for  the  inclination  angle  of  the  twin  jets 
of  the  source  to  the  line  of  sight  of  0  =  66°,  which,  according  to  the  proposed  model,  implies  a  lower  limit 
on  the  central  mass  of  ~  7  x  10 8M©  and  a  lower  limit  for  the  scale  height  of  the  thick  disk  of  ~  12  pc. 
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1.  INTRODUCTION 

Given  that  most  galaxies  harbor  supermassive  black  holes 
at  their  centers  (Richstone  et  al.  1998;  Gebhardt  et  al. 
2000),  and  that  galaxy  mergers  are  common,  binary  black 
holes  should  likewise  be  common.  An  understanding  of  the 
evolution  and  formation  of  these  systems  is  important  for 
an  understanding  of  the  evolution  and  formation  of  galax¬ 
ies  in  general  (Silk  &  Rees  1998;  Merritt  2006).  Theo¬ 
retical  descriptions  of  supermassive  binary  black  hole  sys¬ 
tems  and  their  accretion  disks  have  been  investigated  by 
Hayasaki  et  al.  2007,  2008  and  by  MacFadyen  &  Milosavljevic 
2008. 

Our  ability  to  resolve  the  supermassive  black  holes  in 
any  given  binary  system  depends  on  the  separation  between 
them,  on  their  distance  from  Earth,  and  on  the  resolving 
power  of  the  telescope  used.  It  is  believed  that  the  longest 
timescale  in  the  evolution  of  a  supermassive  binary  black 
hole  system  leading  up  to  coalescence  is  the  stage  in  which 
the  system  is  closely  bound  (~0.1-10  pc),  meaning  that  in 
most  of  these  systems  the  black  hole  pair  can  only  be  re¬ 
solved  by  VLBI  observations  (in  the  case  where  both  black 
holes  are  radio  loud)  which  provide  resolutions  of  milliarc- 
seconds  and  finer.  This  could  explain  why  very  few  such  sys¬ 
tems  have  been  found  (see  review  by  Komossa  2003a  detail¬ 
ing  observational  evidence  for  supermassive  black  hole  bina¬ 
ries). 

Some  source  properties  like  X-shaped  radio  galaxies  and 
double-double  radio  galaxies,  helical  radio  jets,  double-horned 
emission-line  profiles,  and  semiperiodic  variations  in  light 
curves  have  been  taken  as  indirect  evidence  for  compact  binary 

5  Greg  Taylor  and  Ylva  Pihlstrom  are  also  Adjunct  Astronomers  at  the 
National  Radio  Astronomy  Observatory. 


black  holes  though  other  explanations  are  possible.  Some 
wider  systems  have,  however,  been  found  more  directly.  For 
example,  the  ultraluminous  galaxy  NGC  6240,  discovered  by  the 
Chandra  X-ray  observatory ,  was  found  to  have  a  pair  of  active 
supermassive  black  holes  at  its  center  (Komossa  et  al.  2003b), 
separated  by  a  distance  of  1 .4  kpc.  Another  system  that  has  been 
known  for  some  time  is  the  double  AGN  (7  kpc  separation) 
constituting  the  radio  source  3C  75,  which  was  discovered  by 
the  VLA  to  have  two  pairs  of  radio  jets  (Owen  et  al.  1985). 

The  radio  galaxy  0402+379  was  recently  found  to  contain 
two  central,  compact,  flat  spectrum,  variable  components  (des¬ 
ignated  Cl  and  C2),  with  a  projected  separation  of  7.3  pc,  a 
feature  that  had  not  been  observed  in  any  other  compact  source, 
making  this  system  the  most  compact  supermassive  binary  black 
hole  pair  yet  imaged  by  more  than  two  orders  of  magnitude,  with 
an  estimated  system  mass  of  a  few  108  M©  (Maness  et  al.  2004; 
Rodriguez  et  al.  2006). 

Maness  et  al.  (2004)  performed  spectral  line  observations 
of  0402+379  at  1.348  GHz.  Hi  absorption  at  a  redshift  of 
560  km  s_1  from  the  systemic  velocity  of  the  source  (16,489 
±  300  km  s_1;  Xu  et  al.  1994)  was  observed  and  attributed 
to  a  high  velocity  gas  system,  possibly  due  to  a  merger.  Fol¬ 
lowing  this  discovery,  a  spectrum  of  the  H  i  region  was  taken 
in  2004  with  the  Westerbork  telescope  (Morganti  et  al.  2009), 
which  showed  two  components  with  velocities  separated  by 
1000  km  s_1.  In  this  article  we  present  Global  VLBI  observa¬ 
tions  taken  to  improve  the  quality  of  the  Hi  data  on  0402+379, 
with  the  purpose  of  better  understanding  the  origin  of  the  de¬ 
viation  from  the  systemic  velocity  seen  in  the  Hi  gas  in  this 
source. 

At  the  redshift  of  0402+379  of  0.055,  H^  =  15  km  s-1  Mpc-1 , 
and  qo  =  0.5,  a  scale  of  1  mas  =  1.06  pc  is  obtained. 
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Figure  1.  Naturally  weighted  2007  Global  VLBI  image  of  0402+379  at  1.3  GHz.  Contours  are  drawn  beginning  at  3 a  and  increase  by  factors  of  2  thereafter.  In  (b) 
the  synthesized  beam  is  8.15  x  3.72  mas  (shown  in  the  bottom  left  corner).  The  peak  flux  density  is  0.13  Jy  beam-1  and  rms  noise  is  0.02  mJy.  In  (a)  the  image  was 
tapered  and  restored  with  a  circular  25  mas  synthesized  beam  (shown  in  the  bottom  left  corner).  The  peak  flux  density  is  0.46  Jy  beam-1  and  rms  noise  is  0.05  mJy. 


2.  OBSERVATIONS 

2.1.  Global  VLBI  Observations 

Global  VLBI  observations6  were  made  on  2007  March  17  at 
1.3464  GHz.  A  single  intermediate  frequency  with  a  bandwidth 
of  16  MHz  was  observed  with  256  channels  in  both  right  and 
left  circular  polarizations,  resulting  in  a  frequency  resolution  of 
62.5  kHz,  corresponding  to  a  velocity  resolution  of  15  km  s-1. 
Four  level  quantization  was  employed.  The  net  integration  time 
on  0402+379  was  497  minutes. 

Standard  flagging,  amplitude  calibration,  fringe- fitting,  and 
bandpass  calibration  (3C  111  was  used  for  both  gain  and 
bandpass  calibration)  were  followed  in  the  Astronomical  Image 
Processing  System  (AIPS;  van  Moorsel  et  al.  1996).  AIPS 
reduction  scripts  described  in  Ulvestad  et  al.  (2001)  were 
used  for  a  large  part  of  the  reduction.  Spectral  line  Doppler 
corrections  were  also  applied  in  AIPS.  All  manual  editing, 
frequency  averaging  procedures,  imaging,  deconvolution,  and 
self-calibration  were  done  using  Difmap  (Shepherd  et  al.  1995). 
A  clean  cube  was  also  produced  by  first  performing  a  continuum 
subtraction  using  the  task  UVLSF  in  AIPS,  then  using  Difmap 
in  order  to  obtain  a  cleaned  map  for  each  spectral  frequency, 
and  finally  combining  all  the  maps  into  a  cube  with  the  task 
MCUBE  in  AIPS. 

3.  RESULTS 

3.1.  Radio  Continuum 

Figure  1  shows  a  naturally  weighted  1.3  GHz  image  of 
0402+379  from  the  2007  Global  VLBI  observations.  The  source 
consists  of  two  diametrically  opposed  jets  and  a  central  region 
containing  two  active  nuclei  not  resolved  at  this  frequency  (see 
Rodriguez  et  al.  2006  for  a  detailed  multifrequency  study  of 
this  radio  galaxy).  The  northern  jet  is  pointing  in  the  northeast 
direction  whereas  the  southern  jet  is  pointing  in  the  southwest 

6  Global  VLBI  =  Very  Long  Baseline  Array  +  Effelsberg  +  Jodrell  Bank  + 
Westerbork  +  Green  Bank  Telescope  +  Onsala. 


direction.  The  orientation  of  the  source  at  this  frequency  is 
consistent  with  that  seen  by  the  VLA  at  both  1.5  and  5  GHz 
(Maness  et  al.  2004)  and  by  the  VLBA  at  0.3  GHz  (Rodriguez 
et  al.  2006).  The  image  shown  in  panel  (a)  was  tapered  and 
restored  with  a  circular  25  mas  synthesized  beam,  in  order  to 
show  better  the  extended  structure  in  this  radio  galaxy.  The 
source  spans  ~1000  mas  (~1000  pc).  The  image  shown  in  panel 
(b)  was  tapered  and  restored  with  a8. 15x3. 72  mas  synthesized 
beam.  We  see  structure  on  scales  of  ~500  mas  (~500  pc). 

3.2.  H  i  Absorption 

Figure  2  shows  H I  absorption  profiles  taken  from  four  regions 
of  the  source.  The  continuum  has  been  subtracted  from  the 
spectra  using  the  task  UVLSF  in  AIPS,  which  removes  a 
continuum  model  from  the  u ,  v  data  of  all  channels.  The  contours 
are  taken  from  the  2007  Global  VLBI  observations  at  1.3  GHz 
and  the  color  scale  from  the  2005  VLBA  observations  at  5  GHz. 
Two  absorption  lines  are  evident,  which  appear  to  be  at  two 
different  locations  toward  the  southern  jet.  From  this  point  on 
we  will  designate  CW  the  western  component,  where  we  see  the 
stronger  line;  and  CE  the  eastern  component,  where  we  see  the 
weaker  line.  Measurements  of  these  lines  are  given  in  Table  1, 
where  we  also  show  the  peak  opacity  and  the  column  density 
for  both  components,  which  we  calculated  according  to, 

/oo 

z(v)dv(km  s-1) 

-OO 

~  1.8224  x  10187s(K)AD(km  s-1)  ^  r(w),  (1) 

where  Ts  is  the  spin  temperature,  Av  is  the  velocity  resolution, 
and  the  summation  is  over  the  channels  where  we  detect 
absorption  for  CW  and  CE,  respectively. 

We  determined  the  central  velocity  of  the  H I  absorption  lines 
to  be  16,927  d=  7  km  s-1  and  15,856  zb  9  km  s-1  for  CW  and 
CE,  respectively.  From  the  most  recent  reported  redshift  for 
0402+379,  as  measured  from  optical  emission  lines  (Rodriguez 
et  al.  2006),  we  find  that  the  observed  line  for  CW  is  redshifted 
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Figure  2.  Hi  absorption  profiles  taken  from  four  regions  of  the  source  (components  Cl,  C2,  CE,  and  CW).  The  arrow  in  the  center  of  each  spectrum  shows  where  the 
systemic  velocity  is  (Rodriguez  et  al.  2006),  and  the  velocity  resolution  is  15  km  s-1.  The  rms  noise  in  a  single  channel  is  0.14  mJy  beam-1 .  The  contours  are  taken 
from  the  2007  Global  VLBI  observations  at  1.3  GHz  and  are  set  at  1  Oct,  increasing  by  a  factor  of  2  thereafter;  and  the  color  scale  image  is  taken  from  the  2005  VLB  A 
observations  at  5  GHz. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Table  1 

Gaussian  Functions  Fitted  to  Hi  Absorption  Profiles  in  the  Southern  Jet 


Component*1 

Amplitude  (mJy) 

Central  Velocity  (km  s  1 ) 

FWHM  (km  s-1) 

rpeak 

VH i  (cm  2) 

CW 

2.8  ±0.1 

16,927  ±  7 

300  ±  20 

0.025  ±  0.001 

(1.303  ±0.006)  x  1021b 
(7.82  ±0.03)  xl022c 

CE 

1.5  ±0.2 

15,856  ±  9 

170  ±  20 

0.018  ±0.002 

(4.1  ±0.1)  x  1020b 
(2.44  ±0.07)  xl022c 

Notes. 

a  CW  and  CE  refer  to  the  western  and  eastern  jet  components  where  we  find  absorption  lines  respectively. 


b  Assuming  a  spin  temperature  of  100  K. 
c  Assuming  a  spin  temperature  of  6000  K. 

by  370  zb  10  km  s-1  and  for  CE  is  blueshifted  by  700  ± 
10  km  s-1  from  the  systemic  velocity  of  the  source  (16,558 
d=  3  km  s_1).  No  absorption  was  found  at  either  Cl  or  C2  (see 
Figure  2).  We  calculate  a  limit  on  the  H  i  opacity  at  the  location 
of  Cl  and  C2,  and  obtain  less  than  0.17  and  less  than  0.03, 


repectively.  Figure  3  shows  a  map  of  both  the  central  velocity 
and  width  of  the  H  i  absorption  profiles  over  the  source. 

Figure  4  shows  a  velocity  slice  of  the  continuum- subtracted 
cube  (right  panel),  accompanied  by  the  H  i  opacity  distribution 
over  the  source  (left  panel).  This  result  shows  that  either  the  two 
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Figure  3.  (a)  Map  of  the  central  velocity  and  (b)  width  of  the  Hi  absorption  profiles  in  the  source.  Both  maps  were  generated  fitting  Gaussian  functions  at  each  pixel 
where  at  least  a  3 cr  detection  of  a  line  was  found.  The  color  scale  in  both  images  is  in  units  of  km  s_1 .  The  contours  are  taken  from  the  2007  Global  VLBI  observations 
at  1.3  GHz  and  are  set  at  3cr,  increasing  by  a  factor  of  2  thereafter. 


Figure  4.  Velocity  slice  of  the  continuum-subtracted  cube  described  in  Section  3.2  (right  panel),  accompanied  by  the  integrated  H  i  opacity  distribution  over  the  source 
(left  panel),  in  which  we  indicate  where  the  slice  was  taken.  The  map  was  generated  by  combining  a  continuum  image  of  the  source  with  the  continuum-subtracted 
cube.  To  generate  the  opacity  distribution  shown  we  did  a  one-dimensional  zeroth-moment  fitting  to  each  row,  requiring  a  continuum  emission  of  at  least  0.12  mJy 
and  a  line  emission  of  at  least  0.75  mJy  (corresponding  to  ~5<x  respectively).  The  contours  are  taken  from  the  2007  Global  VLBI  observations  at  1.3  GHz  and  are  set 
at  3 cr,  increasing  by  a  factor  of  2  thereafter. 


locations  where  we  find  absorption  are  localized  regions,  rather 
than  being  part  of  a  more  extended  and  perhaps  homogeneous 
structure,  or  our  sensitivity  prevented  us  from  detecting  Hi 
absorption  over  a  broader  region.  In  order  to  explore  this 
question  we  calculated  the  lower  limit  on  the  peak  opacity  across 


the  source  assuming  an  intensity  for  the  line  of  3a ,  where  a  is  the 
rms  noise  in  a  single  channel.  The  result  is  shown  in  Figure  5.  We 
see  that  outside  the  region  where  we  detect  absorption  (compare 
with  Figure  4)  we  would  need  a  peak  opacity  of  at  least  ~0.02, 
comparable  to  the  peak  opacity  found  for  components  CE  and 
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Figure  5.  Calculated  lower  limit  on  the  peak  opacity  across  the  source  assuming 
an  intensity  for  the  line  of  3 o .  The  contours  are  taken  from  the  2007  Global 
VLBI  observations  at  1.3  GHz  and  are  set  at  3 a,  increasing  by  a  factor  of  2 
thereafter. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

CW,  in  order  to  have  a  3a  detection.  The  fact  that  we  did  not 
detect  this  means  that  the  peak  opacity  in  this  region  is  even 
smaller  than  the  lower  limit  found,  suggesting  that  what  we  are 
observing  are  two  different,  localized  clumps,  rather  than  a  more 
extended  and  homogeneous  structure.  The  fact  that  we  do  not 
see  absorption  against  C2  supports  this  scenario. 

We  can  estimate  the  Hi  mass  of  the  clumps  from  the 
calculated  column  densities  (Table  1)  and  assuming  a  spherical 
shape,  according  to  the  following  relation, 

471  9 

m  =  —r  NHlmH ,  (2) 

where  m#  is  the  hydrogen  atom  mass  and  r  is  the  radius  of  the 
clump  estimated  from  the  opacity  distribution  map  shown  in 
Figure  4.  For  CW  we  used  r  ~  12  pc  and  for  CE  r  ~  8  pc.  We 
find  a  mass  ranging  from  ~(7-400)  x  103  M©  for  component 
CW,  and  ~(l-60)xl03  M©  for  component  CE,  assuming  spin 
temperatures  ranging  from  100  K  to  6000  K.  These  masses  are 
bigger  than  typical  values  of  H  i  cloud  masses  of  60  M©  found 
in  the  Milky  Way  (Stil  et  al.  2006). 

4.  DISCUSSION 

The  results  found  from  our  2007  Global  VLBI  observations 
at  1.3  GHz  show  two  Hi  absorption  lines,  one  blueshifted  and 
the  other  redshifted  with  respect  to  the  systemic  velocity.  The 
redshifted  line  (component  CW)  shows  a  FWHM  of  300  zb 
20  km  s_1 ,  and  the  blueshifted  line  (component  CE)  a  FWHM  of 
170  =b  20  km  s_1 .  This  absorption  must  be  from  neutral  hydrogen 
along  the  line-of-sight  to  the  central  continuum  source,  but  this 
gas  could  be  close  to  the  nucleus  or  far  away;  could  be  falling 
in,  outflowing,  or  in  rotation. 


Assuming  a  virialized  cloud  of  ~  104  M©  and  a  radius  of  5 
pc  we  get,  from  the  virial  theorem,  a  velocity  for  the  H  atoms  of 
~  3  km  s_1.  On  the  other  hand,  if  we  use  the  high-end 
temperature  of  104  K,  the  thermal  speed  of  an  H  atom  is  ~  15  km 
s_1.  Thus,  under  these  assumptions,  at  most  ~20  km  s-1  of  the 
line  widths  observed  for  both  CE  and  CW  could  be  accounted 
for,  leaving  a  considerable  balance  which  could  be  attributed  to 
the  proximity  of  these  components  to  the  central  parsecs  of  the 
source,  where  the  gravitational  potential  well  is  deepest. 

Another  possibility  is  that  the  broad  line  widths  and  large 
velocities  for  one  or  both  components  are  due  to  a  jet-cloud 
interaction.  Oosterloo  et  al.  (2000)  studied  the  Seyfert  2  galaxy 
IC  5063,  where  a  strong  interaction  between  the  radio  jet  and 
a  molecular  cloud  of  the  interstellar  medium  (ISM)  is  occur¬ 
ring  at  the  position  of  the  western  radio  lobe  of  the  source. 
The  most  prominent  absorption  feature  found  is  blueshifted 
over  600  km  s-1  with  respect  to  the  systemic  velocity,  a  value 
that  lies  outside  the  range  allowed  by  rotational  kinematics 
of  the  large-scale  Hi  disk.  As  an  effect  of  the  outflow  pro¬ 
duced  by  the  interaction  as  well  as  the  geometry  of  the  source, 
there  will  be  components  moving  both  away  and  toward  the 
observer.  Since  only  absorption  in  those  components  in  front 
of  the  radio  continuum  can  be  observed,  the  blueshifted  com¬ 
ponent  was  the  only  one  detected.  Returning  to  our  case,  it 
is  very  unlikely  that  the  motion  found  in  both  components, 
CW  (redshifted  component)  and  CE  (blueshifted  component), 
which  are  only  a  few  parsecs  apart,  is  due  to  jet-cloud  in¬ 
teractions,  since  this  would  require  a  dramatic  change  in  the 
way  the  southern  jet  is  moving,  on  scales  of  only  a  few  par¬ 
secs. 

We  consider  that  both  CW  and  CE  are  part  of  the  same 
rotating  disk  structure,  of  which  we  only  see  a  couple  of 
clumps,  and  for  the  remainder  of  the  discussion,  we  assume 
that  this  is  the  case  and  propose  a  model  that  reproduces  the 
observations.  We  also  assume  that  the  mass  of  component  Cl 
is  significantly  smaller  than  that  of  component  C2,  thus  not 
affecting  the  stability  of  the  disk.  This  type  of  binary,  where 
the  two  black  holes  have  different  masses,  are  thought  to  occur 
more  frequently,  since  minor  galactic  mergers  are  more  common 
than  major  mergers  in  hierarchical  models  of  galaxy  formation 
(Armitage  &  Natarajan  2002).  We  realize  that  in  a  thin  disk 
the  binary  would  likely  clear  out  all  material  within  twice  the 
radius  of  the  binary  (MacFadyen  &  Milosavljevic  2008),  but 
perhaps  this  process  is  less  efficient  for  a  thick  disk,  or  material 
is  resupplied. 

4.1.  Rotating  Disk 

If  we  assume  that  components  CW  and  CE  are  rotating  about 
C2  in  circular  Keplerian  orbits  lying  on  the  same  plane,  the 
only  way  of  reproducing  the  measured  velocities  (according  to 
our  proposed  model,  explained  below)  is  by  requiring  a  very 
high  central  mass,  of  at  least  ~1010  AT©,  in  disagreement  with 
previous  estimates  (Rodriguez  et  al.  2006)  as  well  as  typical 
supermassive  black  hole  masses. 

Thus,  we  propose  a  model  consisting  of  a  geometrically  thick 
disk,  rather  than  a  geometrically  thin  disk,  rotating  in  a  Keplerian 
circular  orbit,  which  can  be  part  of  a  larger  scale  torus  that 
can  reproduce  the  velocities  which  we  observe  from  the  Hi 
absorption,  at  the  measured  locations.  We  assume  that  CW  and 
CE  are  rotating  in  parallel  planes,  both  offset  from  the  plane 
in  which  C2  lies.  The  distances  between  these  planes  are  free 
parameters  in  our  model.  We  call  dl  the  separation  between 
C2  and  the  center  of  the  circular  orbit  1,  that  CE  is  following; 
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Figure  6.  Configuration  of  the  model  we  are  proposing.  Shown  is  component 
C2,  chosen  as  the  origin,  as  well  as  components  CE,  with  coordinates  (xl, 
_yl),  and  CW,  with  coordinates  (x2,  y2).  The  two  circular  orbits  drawn  are  the 
trajectories  CE  and  CW  would  follow,  respectively,  on  the  rotating  thick  disk, 
shown  in  yellow.  An  inclination  angle  of  ~75°  between  the  rotation  axis  and 
the  line  of  sight  was  used. 


and  d2  the  separation  between  C2  and  the  center  of  the  circular 
orbit  2,  that  CW  is  following.  Figure  6  is  a  cartoon  showing  the 
configuration  of  the  model  we  are  proposing. 

We  assume  that  the  rotation  axis  crosses  component  C2, 
which  is  at  the  systemic  velocity  of  16,558  =b  3  km  s-1 
(Rodriguez  et  al.  2006),  a  natural  assumption  considering  that 
the  jets  of  this  source  seem  to  be  emerging  from  C2,  not  from  C 1 . 
Also,  we  suggest  that  the  two  absorption  lines  that  are  observed 
could  be  explained  by  the  presence  of  two  mass  clumps,  both 
part  of  a  clumpy  geometrically  thick  disk  structure  (an  idea 
that  is  supported  by  the  opacity  distribution  map  obtained, 
Figure  4). 

We  let  d  1  and  d2  vary  over  a  range  of  values  and  calculate, 
for  each  different  combination  (d  \ ,  J2),  what  the  orientation  in 
space  of  the  thick  disk  should  be  in  order  to  best  reproduce 
the  observed  velocities,  via  a  minimization  routine.  After 
obtaining  the  inclination  of  the  thick  disk  for  each  pair  {d  1, 
J2),  we  again  apply  a  minimization  routine  in  order  to  find 
the  combination  (d  1,  d2)  that  best  reproduces  the  measured 
velocities. 

This  same  procedure  was  applied  with  different  values  of 
the  central  mass,  which  allowed  us  to  set  a  lower  limit  on 
this  quantity  of  ~  1.4  x  108  M©  (by  requiring  that  the 
difference  between  the  velocities  predicted  by  the  model  and 
the  observed  velocities  was  smaller  than  la).  See  the  Appendix 
for  more  details  on  the  model.  There  is  no  upper  limit  on 
the  mass  that  we  can  find  using  this  model,  since  for  greater 
masses  we  can  always  find  an  inclination  that  reproduces  very 
accurately  the  measured  values  (the  inclination  angle,  with 
respect  to  the  line  of  sight,  decreases  as  we  increase  the  central 
mass). 

Rodriguez  et  al.  (2006)  made  an  estimate  of  the  central  mass 
of  0402+379,  based  on  an  optical  spectrum  of  the  core  of  this 
radio  galaxy  taken  with  the  Hobby-Eberly  Telescope,  which 
shows  a  red  shoulder  suggesting  two  components  with  a  velocity 
separation  of  ~300  km  s-1 .  At  the  observed  projected  separation 
between  the  nuclei  of  7.3  pc,  an  orbital  velocity  of  300  km  s-1 
implies  a  system  mass  of  1 .5  x  108  M© .  For  this  particular  value 


of  the  central  mass,  our  model  gives  an  inclination  of  the  rotation 
axis  of -85°  with  respect  to  the  line  of  sight,  thus  requiring  the 
source  to  be  very  close  to  the  plane  of  the  sky.  A  value  for  d2 
of  ~  12  pc  is  obtained,  which  sets  a  lower  limit  on  the  scale 
height  of  the  thick  disk.  This  value  is  consistent  with  observed 
parameters  of  supposed  tori  (see  for  example  Peck  et  al.  1999, 
where  a  20  pc  scale  height  torus  is  estimated  for  the  radio  source 
1946+708). 

Considering  the  size  of  our  synthesized  beam  and  the  exten¬ 
sion  of  the  two  regions  where  absorption  is  observed,  we  were 
not  able  to  measure  a  significant  velocity  gradient  within  CW  or 
CE.  The  small  variations  in  velocity  found  are  not  reproduced 
by  our  model.  However,  other  factors  could  influence  these  mass 
concentrations  that  we  are  not  taken  into  account  in  our  model, 
such  as  proper  motions  associated  with  each  clump. 

4.2.  More  Constraints  on  the  Inclination  of  the  Disk 

Assuming  that  C2  is  the  origin  of  the  radio  emission  on 
parsec  scales,  we  can  constrain  the  orientation  of  0402+379. 
In  the  simple  beaming  model  for  simultaneously  ejected  jet 
components  moving  in  opposite  directions,  the  arm  length  ratio 
D  depends  on  the  intrinsic  speed  f  =  v/c  and  the  angle  of  the 
twin  jets  to  the  line  of  sight  0  (Taylor  &  Vermeulen  1997).  The 
arm  length  ratio,  D,  is  given  by 


dN  /  1  +  f  cos  0  \ 

ds  \  1  —  P  cos  0  )  ’ 


(3) 


where  the  apparent  projected  distances  from  C2  (assumed  to 
be  the  origin  of  radio  emission)  are  d n  for  the  northern  jet 
(approaching  side)  and  ds  for  the  southern  jet  (receding  side). 
From  this  we  get  /3cos0  =  0.4  (Rodriguez  et  al.  2006),  a 
result  that  implies  that  the  intrinsic  velocity  must  be  at  least 
0.4c  and  0  must  be  less  than  66°.  Based  on  the  observed 
morphology  of  the  source,  we  can  assume  a  minimum  value  for 
0  of  ~  30°.  According  to  our  proposed  model  (see  Section  4.1) 
these  constraints  on  the  inclination  set  a  lower  limit  on  the 
central  mass  of  the  system  of  ~  7  x  108  Mq,  for  0  =  66°. 
Angles  as  small  as  ~30°  are  allowed,  however,  they  require  a 
very  large  mass  (~1010  M©). 

Figure  7  is  a  cartoon  showing  the  results  obtained  by  our 
model  for  a  central  mass  of  7  x  108M©  ( 0  =  66°),  overlaid  with 
5  GHz  contours  of  0402+379  from  the  2005  VFBA  observa¬ 
tions.  We  can  see  that  the  rotation  axis  aligns  very  well,  within 
~10°,  with  the  jet  axis,  which  was  not  a  constraint  imposed 
by  the  model.  If  we  assume  that  the  center  of  the  orbit  that  Cl 
follows  is  located  at  the  position  of  C2,  and  also  that  component 
C2  is  significantly  more  massive  than  Cl,  we  obtain  an  orbital 
velocity  for  Cl  of  ~  580  km  s-1  (for  a  mass  of  7  x  108  M©  and 
0  =  66°).  This  is  consistent  with  the  current  limit  from  VFBI 
monitoring  of  less  than  26,400  km  s-1 . 


5.  CONCFUSION 

Global  VFBI  observations  at  1.3  GHz  were  performed  on  the 
radio  galaxy  0402+379,  the  most  compact  supermassive  binary 
black  hole  pair  yet  imaged.  Two  absorption  lines  were  found 
toward  the  southern  jet  of  the  source,  one  redshifted  by  370  zb 
10  km  s-1  and  the  other  blueshifted  by  700  ±10  km  s-1  with 
respect  to  the  systemic  velocity  of  the  source. 

A  model  consisting  of  a  thick  disk,  of  which  we  only  see 
a  couple  of  clumps,  was  developed  in  order  to  reproduce  the 
velocities  measured  from  the  Hi  absorption  profiles.  These 
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Figure  7.  Results  obtained  from  the  rotating  thick  disk  model,  projected  on  the 
x-y  plane  (R.A.  and  decl.  plane),  overlaid  with  5  GHz  contours  from  the  2005 
VLB  A  observations.  The  blue  section  of  the  rotation  axis  is  pointing  towards  us, 
whereas  the  red  section  is  pointing  away.  For  a  mass  of  7  x  108M©  an  inclination 
angle  of  ~66°  between  the  rotation  axis  and  the  line  of  sight  is  obtained.  The 
position  of  component  C2,  chosen  as  the  origin,  is  shown,  as  well  as  the  position 
of  Cl,  CE,  and  CW.  The  two  green  circular  orbits  drawn  are  the  trajectories 
CE  and  CW  would  follow,  respectively,  on  the  rotating  thick  disk.  The  yellow 
circular  orbit  is  the  trajectory  Cl  would  follow  assuming  C2  as  the  center  of  the 
orbit. 


clumps,  components  CW  and  CE,  rotate  in  circular  Keplerian 
orbits  around  an  axis  that  crosses  one  of  the  supermassive  black 
holes  of  the  binary  system  in  0402+379,  component  C2.  We 
found  an  upper  limit  for  the  inclination  angle  of  the  twin  jets  to 
the  line  of  sight  of  0  =  66°,  which,  according  to  the  proposed 
model,  implies  a  lower  limit  on  the  central  mass  of  -7x10  8M0 
and  a  lower  limit  for  the  scale  height  of  the  thick  disk  of 
~  12  pc. 

Facilities:  VLB  A,  Effelsberg,  Jodrell  Bank,  Westerbork, 
Green  Bank  Telescope,  Onsala  25  m  telescope 


is  unitary  we  have, 

n1x+n2y+n2z  =  1.  (Al) 

According  to  Kepler’s  Laws,  a  body  moving  around  a  mass 
M,  at  a  radius  R ,  has  a  velocity  given  by 


Thus,  the  angular  velocity  is  given  by, 

co  =  VGM(R)~3/2.  (A3) 

If  ifr  is  the  angular  velocity  vector  of  our  rotating  thick  disk, 
and  is  any  point  on  the  disk,  then  the  tangential  velocity 
is  given  by 

V  =  7$  x  7t  (A4) 

where  7$  =  am  and  ~K  =  T*  —  ~ct ,  where  is  the  separation 
between  the  origin  and  the  center  of  the  circular  orbit  we  are 
considering,  and  is  the  position  vector  of  the  point  we  are 
considering  (in  our  case,  the  coordinates  of  either  CE  or  CW). 
Thus,  we  have, 

=  (RX1  Ry ,  Rz )  =  (jc  +  dnx ,  y  +  driy ,  z  +  dnz).  (A5) 

Let  the  plane  x-y  be  the  plane  of  the  sky  and  z  be  the 
line-of- sight  direction.  We  want  to  know  the  z-component 
of  the  tangential  velocity  ~V^,  because  that  component  is  the 
one  we  can  measure  (from  our  absorption  lines),  which  is 
given  by 

vz  =  V  -z  =  co(nxRy  -  nyRx).  (A6) 

Finally,  we  can  also  use  the  fact  that  7%  and  7?  have  to  be 
orthogonal,  that  is, 

=  0  =>-  nxx  +  nyy  +  nzz  +  d  =  0.  (A7) 

We  can  now  put  everything  together  and  get  an  expression 
for  Vz, 

Vz  =  V GM  [nx(y  +  dny )  —  ny(x  +  dnx)\ 

X  [(x  +  dnx)2+(y  +  dny)2+(z  +  dnzf]~V4 ,  (A8) 

where 
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APPENDIX 

ROTATING  DISK  MODEL 

Figure  6  illustrates  the  configuration  of  the  model  we  are 
proposing.  We  show  component  C2,  which  is  chosen  as  the 
origin,  as  well  as  components  CE,  with  coordinates  (xl,  yl),  and 
CW,  with  coordinates  (x2,  y2).  The  two  circular  orbits  drawn  are 
the  trajectories  CE  and  CW  would  follow,  respectively,  on  the 
rotating  thick  disk.  Also,  we  show  the  projection  of  the  system 
on  the  plane  of  the  sky,  since  that  is  what  we  actually  observe. 

Let  us  assume  our  disk  has  an  arbitrary  orientation  in  space, 
specified  by  the  normal  unitary  vector  n  —  (%,  ny,  nz).  Since  n 


xnx 


z  = 


yny 


and  n7  = 


= 71  -»? 


(A9) 


Thus,  we  have  an  equation  that  gives  us  the  z-component  of 
the  tangential  velocity  of  any  point  with  a  circular  orbit  in  our 
rotating  disk,  in  terms  of  the  v-  and  y-components  of  the  normal 
unitary  vector  of  the  disk,  nx  and  ny,  and  in  terms  of  the  v-  and 
y-components  of  any  point  in  the  disk. 

There  are  two  locations  in  our  source  where  we  can  measure 
the  velocity.  Let  us  define  Vz\t0  and  VZ2,0  as  the  observed 
velocities  of  CE  and  CW  respectively;  and  Vz\tC  and  VZ2,C  the 
calculated  velocities  according  to  Equation  (A8).  Finally  let  us 
define  the  function  that  we  will  minimize  in  order  to  obtain 
the  inclination  of  the  thick  disk  that  reproduces  the  measured 
velocities  best, 


|Vzl,o-Vd,cl  ,  \Vz2,o  -Vz2, cl 
\VzU\  +  I  Vz2,o  I 


(A10) 
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For  a  particular  value  of  the  central  mass  of  the  system, 
M,  we  went  through  the  following  routine  in  order  to  find 
the  inclination  of  the  thick  disk,  more  precisely  nx  and  ny 
(since  nz  can  be  determined  in  terms  of  these  two  quan¬ 
tities):  we  calculated  the  function  J  for  every  combination 
of  ( nx ,  ny ),  with  the  constraint  that  n  is  unitary.  By  mini¬ 
mizing  J  we  found  the  best  pair  (, nx ,  ny).  Since  d  1  and  d2 
are  also  free  parameters,  we  let  them  vary  over  a  range, 
so  for  every  combination  (d\ ,  d2)  we  followed  the  proce¬ 
dure  mentioned  above  and  found  the  best  pair  (, nx ,  ny).  Af¬ 
ter  having  done  this  for  every  combination  (d  1 ,  d2),  again 
by  a  minimization  routine,  we  found  the  best  pair  (d  1,  J2). 
Thus,  for  a  particular  value  of  the  central  mass  we  were  able 
to  find  the  inclination  of  the  disk  that  reproduced  the  ob¬ 
served  velocities  best.  This  same  routine  was  also  followed 
for  different  values  of  the  central  mass,  which  let  us  find  a 
lower  limit  on  this  quantity  of  ~  1.4  x  10 8MQ,  by  requir¬ 
ing  that  the  difference  between  the  velocities  predicted  by 
the  model  and  the  observed  velocities  was  smaller  than  la, 
with 


a  = 


j 


(VzUo  -  Vzlcf 

A? 


+ 


(All) 


where  Ai  and  A2  are  the  errors  in  the  measurements  of  Vzii0  and 
VZ2 ,0  respectively. 
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